Parametric excitation of a magnetic nanocontact by a microwave field 
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We demonstrate that magnetic oscillations of a current-biased magnetic nanocontact can be para- 
metrically excited by a microwave field applied at twice the resonant frequency of the oscillation. 
The threshold microwave amplitude for the onset of the oscillation decreases with increasing bias 
current, and vanishes at the transition to the auto-oscillation regime. The parametrically excited 
oscillation mode is the same as the one in the auto-oscillation regime, enabling studies of both the 
passive and the active dynamics of the oscillator. Theoretical analysis shows that measurements of 
parametric excitation provide quantitative information about the relaxation rate, the spin transfer 
efficiency, and the nonlinearity of the nanomagnetic system. 

PACS numbers: 05.45.Xt,76.50.+g,75.75. Jn,85.70.Ec 



Several novel magnetic nanodevice architectures have 
been recently proposed for the applications in informa- 
tion technology JEUll , and for generation 0,0], sensing [U 
and processing [6|,LZ| of electromagnetic signals. Their im- 
plementation critically depends on our ability to quanti- 
tatively characterize and control the dynamical charac- 
teristics of nanomagnets. One of the most significant 
recent developments that provided insight both into the 
dynamical properties and the mechanisms of excitation 
of nanomagnetic systems is the spin-torque ferromagnetic 
resonance technique (ST-FMR) [8l4l0l|. an extension of 
the ferromagnetic resonance technique commonly utilized 
for the characterization of magnetic materials . 

In the ST-FMR method, a microwave current with fre- 
quency f e close to the resonance frequency f of nano- 
magnet is applied to the nanomagnetic device. A dc volt- 
age is produced by mixing of the microwave current with 
the signal generated by the dynamical response of the 
nanomagnet. By modeling the dependence of this volt- 
age on the applied microwave frequency, one can extract 
information about the characteristic frequencies, relax- 
ation rates, and the spin-polarization of electrical current 
in the nanomagnetic system. 

Another method previously developed for the studies 
of magnetic materials is the parametric pumping spec- 
troscopy, which utilizes microwave-frequency modulation 
of the applied field to excite magnetic dynamics llj, [12 1 . 
This technique provides information complementary to 
FMR about the dynamical properties of magnetic materi- 
als. For instance, FMR measurements can be affected by 
simultaneous excitation of several dynamical modes, re- 
sulting in jumps of the resonant frequency and linewidth 
broadening Q. In contrast, parametric excitation has a 
threshold nature, providing information about a single 
excited mode at driving signals that are not too large. 

In this Letter, we report the first observation of para- 
metric excitation of a nanomagnet by a microwave mag- 
netic field, applied at frequency f e equal to twice the 



resonance frequency /o of the nanomagnet. Although 
our nanomagnetic system is a nanocontact on a spatially 
extended magnetic film that has a continuous excitation 
spectrum, we demonstrate that only one dynamic mode 
characterized by the lowest damping is parametrically 
excited, enabling an accurate determination of the spe- 
cific parameters of this mode. We show that the depen- 
dence of parametric excitation on the driving frequency is 
strongly asymmetric, which is caused by the nonlinearity 
of the studied dynamical sysytem. We also demonstrate 
that all the important features of our observations can 
be quantitatively described by the analytical model of a 
non- autonomous nonlinear oscillator T3, HI- 

Based on our observations and the developed theory, 
we propose a simple quantitative method for the char- 
acterization of magnetic nanoelements. We demonstrate 
that by measuring the threshold and frequency range of 
parametric excitation, it is possible to determine such 
important parameters as damping, spin-polarization effi- 
ciency, and coupling coefficient to the microwave signal. 
In addition, by measuring the frequency range of para- 
metric synchronization in the auto-oscillation regime, one 
can independently determine the dynamic nonlinearity of 
the nanomagnet. A significant advantage of the proposed 
parametric approach over the ST-FMR technique is pro- 
vided by the ability to directly measure the induced oscil- 
lation by spectroscopic techniques, without any interfer- 
ence from the pumping signal whose frequency f e w 2/o 
is significantly higher than /o. 

Measurements of parametric excitation were 
performed in nanocontact devices with structure 
Cu(40)Py(3.5)Cu(8)Co 70 Fe 30 (10)Cu(60), fabricated 
on sapphire substrates with electrical leads patterned 
into coplanar microstrip lines. Here, thicknesses are in 
nanometers, and Py = NisoFe2o- The polarizing CoFe 
layer and part of the Cu(8) spacer were patterned into 
an elliptical shape with dimensions of 100 nmx50 nm. 
The free Py(3.5) layer was left extended with lateral 



2 



microstrip H) a 
H h I**- 



polarizer- 
free layer 



10 



I 

3 



o 




% (d) 











\o J 






\° / 




2.0 2.5 


3.0 




I (mA) 






FIG. 1: (Color online) (a) Schematic of the studied point 
contact device with a microstrip generating microwave field 
h e . (b)-(d) generation frequency, power, and linewidth vs 
bias current I, at h e = 0. Symbols are data, curves in (b), 
(c) are theoretical fits, and in (d) - a calculation based on the 
nonlinear auto-oscillator model (Eq. (95) in Ref. [13]), with 
parameters determined from (b) and (c). 



dimensions of several micrometers (Fig. [IJa)), resulting 
in a device geometry similar to the point contacts studied 
before [i, 15|. The bias field H = 1.1 kOc was oriented 
in the device plane, at 45° with respect to the easy axis 
of the nanopatterned CoFe layer. Measurements were 
performed at 5 K. We show data for one of two devices 
tested with similar results. 

The dynamical properties of nancontacts were char- 
acterized by measurement of auto-oscillation induced at 
bias current I > I c in the absence of the external driv- 
ing signal. Here, I c = 2.0 mA is the critical current 
for the onset of auto-oscillation. The dependence of the 
auto-oscillation frequency /o on current exhibited a slight 
increase just above I c , and a decrease at I > 2.7 mA 
(Fig. D^b)). The generated power monotonically in- 
creased with I (Fig.[IJc)), while the linewidth (Fig.[ljd)) 
exhibited a non-monotonic behavior consistent with the 
effects of nonlinarity on thermal line broadening |16l Il7j . 

The pumping microwave field h e 1 1 H was generated by 
a microwave current i m applied to a Cu microstrip fabri- 
cated on top of the nanocontact and electrically isolated 
from it by a SiC>2(50) layer (Fig. Ufa)). The dependence 
of the microwave field on the ac current was caliabrated 
by a procedure described elsewhere [ljl [l8[ . To paramet- 
rically induce oscillations, a microwave field at frequency 
/ e « 2/o was applied to the device. We found that even 
the largest field h e — 35.6 Oe rms in our measurements 
was below the oscillation threshold at I = 0. However, 
we were able to induce oscillations by simultaneously ap- 
plying h e and a subcritical bias current I > 1 mA that 
partially compensated the damping (l9| . 

In spectroscopic measurements performed at differ- 
ent values of / e , the oscillations appreared near f e = 
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FIG. 2: (Color online) (a)-(c) Parametric excitation in the 
subcritical regime, at / = 1.7 mA, and (d) parametric syn- 
chronization in the supercritical regime, at / = 2.3 mA. (a) 
Spectra of the parametrically excited oscillations at the la- 
beled values of f e , at h e = 12.6 Oe rms. Curves are offset 
for clarity, (b) Full width at half maximum (FWHM) of the 
oscillation peaks under the same conditions as in (a) ; (c) Fre- 
quency boundaries of the parametric excitation region vs h e . 
(d) Frequency boundaries of the parametric synchronization 
region vs h e . Symbols in (c) and (d) are data, curves - calcu- 
lations using Eq.Q and Eq.®, respectively, as described in 
the text. 



16.95 GHz 2/o(I c ) (Fig. Ufa)), independently of the 
bias current 1 mA < I < I c = 2 mA or h e . Here, 
fo{Ic) = 8.465 GHz is the oscillation frequency just above 
the auto-oscillation onset, at h e = (Fig. [TJb)). There- 
fore, we can conclude that the parametrically excited os- 
cillation mode is the same as the one generated at the 
onset of auto-oscillation, for a wide range of / < I c and 
h e . This observation enabled us to directly compare the 
dynamical characteristics of the device extracted from 
the parametric excitation to the properties known from 
the measurements of autonomous dynamics (Fig[TJ) . 

The amplitude of the driven oscillation exhibits a max- 
imum near zero detuning A/ = f e /2 — fo(I c ) (Fig.[2ja)), 
while the linewidth has a minimum near this point 
(Fig. E[b)). The frequency of the oscillation is exactly 
equal to / e /2. The oscillation completely vanishes at fre- 
quencies f e < f e ,miri and f e > f e , max - The frequency 
range of the parametric excitation is proportional to the 
driving amplitude h e above a threshold value h t h = 8 Oe 
(dots in Fig. [U[c)). We show below that the slope of this 
dependence is determined mainly by the intrinsic damp- 
ing, in agreement with the general properties of paramet- 
ric excitation. 

In the supercritical regime (7 > I c ), the oscillation was 
observed for all values of f e . At f e w 2/o, it became syn- 
chronized with the microwave field, similarly to the mag- 
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FIG. 3: (Color online) (a) Symbols: measured microwave os- 
cillation power us pumping frequency, at h e = 12.6 Oe rms 
(filled symbols) and h e = 22.5 Oe rms (open symbols) in the 
subcritical regime, at I = 1.7 mA. Curves are calculations 
using Eq. ([2]). (b) Threshold microwave amplitude vs bias 
current for several values of detuning: A/ = (filled symbols 
and solid curve), A/ = 150 MHz (crosses and dashed curve), 
and A/ = —150 MHz (open symbols and dotted curve). Sym- 
bols are data, curves are calculations using Eq. @ for A/ = 0, 
150 MHz, and Eq. © for A/ = -150 MHz. 



netic nanopillars as described in Ref. [18j. The depen- 
dence of the synchronization boundaries on h e (Fig.[^d)) 
appears to be similar to the parametric excitation data 
(Fig. [D^d)). However, we show below that the synchro- 
nization interval is determined by the dynamic nonlinear- 
ity of the device rather than damping. In contrast to the 
parametric excitation, the synchronization is observed at 
any h e , i.e. h th = 0. 

The dependence of the oscillation power P on f e is 
strongly asymmetric with respect to the sign of A/ 
(Fig-ISI a ))- At A/ < 0, the oscillation amplitude grad- 
ually decreases to zero with increasing magnitude of de- 
tuning, while at A/ > it initially increases with de- 
tuning and then abruptly decreases to zero. The depen- 
dence of the excitation threshold h t h on / (Fig. E^b)) 
is also asymmetric with respect to the sign of A/. For 
A/ = 150 MHz, the threshold linearly decreases with /. 
For A/ = -150 MHz, the threshold closely follows the 
A/ = 150 MHz values at / < 1.5 mA, while at larger / 
the decrease becomes slower. 

To understand the origin of these unusual features of 
parametric excitation of the nanomagnetic oscillator, we 
utilize the model of a driven nonlinear oscillator devel- 
oped in Refs. [13, 14 1. The state of the oscillator is 



characterized by a dimensionless complex amplitude c(i), 
whose evolution is determined by 

dc 
dt 



+ iuj(p)c + T(I,p)c = Vh e e 



(1) 



Here, p = |c| 2 is the dimensionless oscillation power, 
u e = 2Trf e , u)(p) = 27r/o(/ c )[l + £,{v)] is the power- 
dependent auto-oscillation frequency, V is the coupling 
to the driving field h e , and T(I,p) is the total effective 
damping given by the difference between the natural pos- 
itive damping T + (p) = To[l + and the negative 



damping r_(/,p) = ro(/// c )(l — p) caused by the spin- 
polarized current I. 

The functions and rj{p) characterize the nonlinear- 
ities of the oscillation frequency and the natural damping, 
respectively. They were determined by fitting the data 
of Figs. Dlb), (c) with Eq. ([T]), with the right-hand side 
taken to zero. The applicability of the model Eq. ((T|) to 
the studied magnetic nanocontacts was verified by an in- 
dependent calculation of the generation linewidth. The 
known functions £(p) and rj{p) allowed us to determine 
the dimensionless power-dependent nonlinearity coeffi- 
cient v{I,p) = [duj(p)]/dp]/[dr(I,p)/dp], and to analyt- 
ically calculate the nanocontact generation linewidth in 
the active regime (Eq. (95) in Ref. [13]). The temper- 
ature for this calculation was elevated to T = 10 K to 
account for the Ohmic heating of the nanocontact by the 
bias current [20|. The calculation shows a good agree- 
ment with the data over a large range of I (curve in 
Fig. Hid)), supporting the applicability of our nonlinear 
oscillator model to the studied magnetic point contacts. 

To analyze the mechanisms of parametric excitation, 
we note that Eq. (|TJ) admits a synchronous solution in the 
form c(t) = */pe _ia,ei/ ' 2+l ^ 1 , where p > is the oscillation 
power, and ip is its phase relative to the driving signal. 
Substituting this expression into Eq. ([1]) and multiplying 
both sides by their complex conjugates, we obtain an 
implicit expression for p 
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+ T z (I,p)=V'h 



(2) 



A solution to this equation exists for h e > hth, where the 
excitation threshold hth corresponds to the minimum of 
the left-hand side of Eq. . Depending on the param- 
eters of the model, this minimum can occur either at 
p = or at some finite power pj > 0. The former case 
corresponds to the "soft" regime of parametric excita- 
tion, with the oscillation power p gradually increasing 
from p = with increasing h e > hth, whereas the lat- 
ter case describes the "hard" regime characterized by an 
abrupt jump of p from zero to a finite value pf. 

Solving Eq. ^ in the "soft" regime, we obtain the 
threshold microwave amplitude 



Vh' th = ^/A^+Tj, (3) 
and the boundaries of the parametric excitation region 



■■h 



2w n ± 2\ V 2 hl 



r 2 



(4) 



where Aw = 2-7rA/ = w e /2 — ujq is the linear frequency 
detuning and Tj = Tq(1 — I/I c ) is the linear damping 
reduced by the subcritical current I. 

The linear oscillation frequency Wo, linear damping 
rate Tq, parametric coupling coefficient V, and critical 
current I c can be determined using Eq. Q from the de- 
pendence of the parametric excitation frequency interval 



4 



on the driving field. From the data of Fig. [2jc), we ob- 
tained w /2tt = f = 8.475 GHz, r = 1.31 ns" 1 (corre- 
sponding to the Gilbert damping coefficient a G = 0.015), 
V = 2tt • 3.31 MHz/Oe, and J c = 1.99 mA. 

The "hard" regime of parametric excitation occurs 
when the power-dependent detuning Aw(p) = w e /2 — 
lo(jp) decreases with the increase of the oscillation power 
p. Thus, for a given sign of the nonlinearity coefficient v, 
the "hard" excitation takes place only on one side of the 
resonance uj e /2 = ujq. The threshold h" b in the "hard" 
regime is approximately given by 



\Auj + is(I, O^jI 



(5) 



The experimental dependence of the auto-oscillation 
frequency on the bias current exhibits an initial increase 
at I c < I < 2.7 mA, thus v > 0, and therefore the 
"hard" regime of parametric excitation only at Auj > 0. 
The resulting asymmetry of the dependence of the os- 
cillation power on detuning is clearly seen in Fig. |3j[a). 
The difference between the "soft" and the "hard" regimes 
is also illustrated in Fig. [31(b): in the "hard" regime at 
A/ = 150 MHz, the parametric threshold linearly de- 
creases with the bias current (Eq. ((5$), while in the "soft" 
regime at A/ = — 150 MHz it exhibits a nonlinear depen- 
dence on / (Eq. 0). 

At / > I c , the parametric excitation is replaced by 
the parametric synchronization characterized by the syn- 



chronization index r — uj e /uj = 2 
expression for the frequency range of the parametric syn- 
chronization can be obtained by using Taylor expansions 
of w(p) and T(I,p) around the free-running power p = pa 
in Eq. J3]), yielding 



<i, max / min 



2co(p ) ± 2y/l + v 2 (I,p Q )Vh e 



(6) 



This expression reasonably well describes the experimen- 
tal data in Fig. H][d) , proving that measurements of the 
parametric synchronization range can be used for the in- 
dependent determination of the nonlinear coefficient v in 
nanomagnetic oscillators. This result is especially im- 
portant for the studies of other potentially more compli- 
cated dynamical systems such as magnetic nanopillars, 
where in some cases several magnetic modes can be ex- 
cited simultaneously 16j, affecting the measurements of 
nonlinearity in the autonomous regime. 

To summarize, we have reported the first observation 
of asymmetric parametric resonance in a current-biased 
magnetic nanocontact, and demonstrated that this phe- 
nomenon can be utilized to determine dynamical prop- 
erties of magnetic nanoelements. We have also demon- 
strated that the general model of a nonlinear oscillator 
[l3| provides a quantitative description of the observed 
autonomous and non-autonomous (driven) nanomagnet 



dynamics. The parametric measurements can be utilized 
as an efficient characterization technique complementary 
to the ST-FMR method for the studies of the dynamical 
properties of nanoscale magnetic systems. 
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